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ABSTRACT 

Utilizing the Two Micron All Sky Survey (2MASS) Second Incremental Data Release 
Catalog, we have retrieved near-IR magnitudes for several hundred hot subdwarfs (sdO 
and sdB stars) drawn from the Catalogue of Spectroscopically Identified Hot Subdwarfs 
(Kilkenny, Heber, & Drilling 1988, 1992). This sample size greatly exceeds that of pre- 
vious studies of hot subdwarfs. Examining 2MASS photometry alone or in combination 
with visual photometry (Johnson B For Strdmgren uvby) available in the literature, we 
show that it is possible to identify hot subdwarf stars that exhibit atypically red IR 
colors that can be attributed to the presence of an unresolved late type companion. 
Utilizing this large sample, we attempt for the first time to define an approximately 
volume limited sample of hot subdwarfs. We discuss the considerations, biases, and 
difficulties in defining such a sample. 

We find that, of the hot subdwarfs in Kilkenny et al., about 40% in a magnitude 
limited sample have colors that are consistent with the presence of an unresolved late 
type companion. Binary stars are over-represented in a magnitude limited sample. In an 
approximately volume limited sample the fraction of composite-color binaries is about 
30%. 

Subject headings: binaries: general — stars: early-type — stars: horizontal-branch - 
subdwarfs — infrared: stars 

1. Introduction 

Hot subdwarf stars as a class have effective temperatures exceeding ~20,000 K, with surface 
gravities higher and luminosities lower than main sequence stars of the same temperature. Spectro- 
scopically, hot subdwarfs are divided into two major subtypes: subdwarf B (sdB) and subdwarf O 
(sdO). SdB stars have hydrogen-dominated atmospheres with T e g « 24, 000-30,000 K and exhibit a 
relatively small dispersion in luminosity (Heber 1986; Saffer et al. 1994; de Boer et al. 1997). They 
often exhibit weak lines of Helium, Carbon, and other atomic species in comparison with Popula- 
tion I main sequence stars. SdO stars have helium-enriched atmospheres with T e g- > 30, 000 K and 
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exhibit a wider range of log g and luminosity than the sdB stars (Heber 1986; Saffer et al. 1994; 
de Boer et al. 1997). Hot subdwarf stars (and sdB stars in particular) are the primary source of 
the "UV upturn" (or "UVX") seen in elliptical galaxies and spiral galaxy bulges (e.g., Greggio & 
Renzini 1990; Brown et al. 1997, 2000; O'Connell 1999). 

Based on their temperatures and magnitudes, Humason & Zwicky (1947) proposed that the 
sdB stars are the field equivalent of Baade's (1944) population n B stars (now known as "extended 
horizontal branch 1 ," or "EHB," stars) which are seen in some globular cluster HR diagrams at the 
blue end of the horizontal branch. As such, sdB stars in the field are understood to be core helium 
burning objects with very low hydrogen envelope masses (M env < 0.02 M Q ) and total masses of 
M 0.5 — 0.55 M Q (Caloi 1972; Saffer et al. 1994). Some sdO stars have colors and luminosities 
that are consistent with the hottest part of the EHB (least massive H envelopes). However, there 
are a number of potential evolutionary paths, e.g. post-AGB evolution, that lead through the rather 
broad region of the HR diagram occupied by sdO stars. All hot subdwarfs, regardless of their prior 
evolution, are thought to be direct progenitors of white dwarfs, although they constitute only a 
small fraction of all stars that evolve to become white dwarfs. 

While the current evolutionary status of sdB stars is relatively well understood, how they 
actually arrived at the EHB with such low envelope masses is not understood. The spectroscopically 
inferred masses of field sdB stars have a very small dispersion. If they evolved directly from red 
giant progenitors, then their initial masses are constrained to be < 2 M , leading to a requirement 
that they lose up to ~1.5 M Q on the red giant branch before or right at core helium ignition 
(Saffer et al. 1994). One scenario for such extreme mass loss suggests that sdB stars lost significant 
envelope mass due to binary interactions on the red giant branch (Mengel, Norris, & Gross 1976). 
If this were the sole channel for forming sdB stars, all sdB stars would be expected to be found in 
close binary systems — that is systems where the current separation is of the same order as the 
radius of the sdB progenitor star near the tip of the red giant branch, or R ~ 100 R0 ~ 0.5 AU 
(Maxted et al. 2001). Companions to sdB stars have been reported by numerous surveys, involving 
a variety of photometric and spectroscopic techniques, although these often have been of limited 
scope. These include: searches for flux excesses at long wavelengths, or flat energy distributions 
(e.g., Allard et al. 1987; Allard et al. 1994; Thejll, Ulla, & MacDonald 1995; Ulla & Thejll 1998); 
radial velocity variations (e.g., Maxted et al. 2001; Saffer, Green, & Bowers 2001); the presence of 
a G-band, Ca II absorption, or other spectral features indicative of late type stars (e.g., Ferguson, 
Green, k Liebert 1984; Bixler, Bowyer, & Laget 1991; Theissen et al. 1995; Jeffery & Pollacco 1998; 
Saffer, Green, & Bowers 2001; Maxted et al. 2001); and even direct high-resolution imaging (e.g., 
Thejll, Theissen, & Jimenez 1994; Heber et al. 2002). 

Some of these survey reports have included an attempt to deduce the fraction of sdBs that are 
in binaries. The results have ranged from 50 to 100%, depending on differing specific assumptions 
made about the nature and physical parameters of the sdB stars and/or their companions. For 



lr The term "extended horizontal branch" was coined by Greenstein & Sargent (1974). 
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example, the evolutionary status of the companion stars is typically not known. Some authors 
have preferred main sequence companions and less luminous sdB stars (i.e., Allard et al. 1994; 
Ferguson et al. 1984; Bixler et al. 1991), whereas others have preferred slightly evolved, sub-giant 
or giant companions and more luminous sdB stars (i.e., Jeffery & Pollacco 1998; Theissen et al. 
1995). The correction for "unseen" fainter companions (such as dM stars) in determining the total 
binary fraction is greatly affected by this choice. Thus, if the observed (G-K type) companions 
in composite systems are subgiants, the number of fainter unobserved companions (extrapolated 
using an assumed initial mass function) must be very large, and the resulting binary fraction is 
close to 100%. On the other hand, if the G-K companions are main sequence stars, the number of 
fainter unobserved companions is much smaller (because the lower inferred luminosity for the hot 
subdwarf is less effective at masking the light of such companions). The result is a binary fraction 
as low as 50%. However, only the radial velocity or imaging studies give information about whether 
the sdBs that do have companions are in close binary systems. 

The systematic survey of the entire sky that was carried out by the Two Micron All Sky 
Survey 2 (2MASS) provides the opportunity to further explore the question of whether all sdB stars 
(and sdO stars, to the extent that they represent the EHB) are in binary systems. This can be done 
by combining cataloged photometry at visible wavelengths with the 2MASS near-infrared imaging 
photometry. In addition to comparing visible and near-IR magnitudes, near-infrared colors for hot 
subdwarfs (e.g., J — K$) can be used on their own as a second indicator of a binary system. In 
this paper we collect readily available visible and near-IR flux measurements of hot subdwarfs and 
identify those whose colors indicate the presence of a late type companion. We thus determine 
the fraction of subdwarfs that exist in composite spectrum binaries. (Note that not all classes 
of companion star will be identified this way, in particular white dwarf companions.) We extract 
visible photometry from an existing catalog of hot subdwarf stars, and we use IR photometry 
obtained from the 2MASS Second Incremental Data Release (2IDR) Catalog 3 . We look for IR 
colors that are "redder" than what would be expected for a normal hot subdwarf and which can 
be attributed to the presence of an unresolved late type companion. This method is similar to that 
used by Thejll et al. (1995) and Ulla Sz Thejll (1998), but our sample is much larger and the IR 
photometry is more precise. 

In §2 we describe our sample of hot subdwarfs including: improved coordinates, retrieval of 
2MASS infrared photometry, availability and sources of visual photometry, and adopted classifica- 
tions for objects in the initial sample. In §3 we examine and model a magnitude limited sample of 
hot subdwarfs using two color diagrams and the distribution in J—Ks color. In §4 we introduce an 
approximately "volume limited" sample of hot subdwarfs and compare it to the magnitude limited 
sample. Finally, in §5 we summarize our findings. 



ht t p : / /www .ipac.caltech.edu/2mass/ 
3 http: / /www. ipac.caltech.edu/2mass /releases/second /index. html 
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2. The Sample 

For this investigation we studied hot subdwarf stars listed in the Catalogue of Spectroscopi- 
cally Identified Hot Subdwarfs (Kilkenny, Heber, & Drilling 1988) as updated and expanded in an 
electronic version by Kilkenny, c. 1992, kindly made available to us by D. Kilkenny. This catalog 
(hereinafter referred to as KHD) contains 1527 entries, of which 302 entries were added since the 
1988 publication. Hot subdwarfs from the PG survey of UV-excess objects (Green, Schmidt, & 
Liebert 1986) are an important component of the catalog (1017 objects with PG designations, either 
discovered or recovered by that survey). More recent or ongoing surveys, including the Hamburg 
Schmidt (HS) and Edinburgh-Cape (EC) surveys are not included. In the catalog sub-sample of 
593 objects that is the subject of this paper, about 84% of the stars are to be found in the 1988 
printed version of KHD, and about half of the sub-sample stars added to KHD between the 1988 
and 1992 versions are from the Carnegie Southern Observatory H&K survey (Beers et al. 1992). 
About 71% of the entries in KHD have northern Declinations, the same percentage as in the sub- 
sample studied here. The sub-sample is thus broadly representative of the parent catalog, for stars 
bright enough to be included in the 2MASS point source database. A number of important selection 
effects bias the parent catalog's contents, however, since it is compiled from many different sources 
each with its own criteria for inclusion or discovery (see also §3.3), but KHD is the most complete 
single compilation of field hot subdwarfs available, especially for the northern sky and for relatively 
bright objects. 

While the KHD catalog contains all varieties of hot subdwarfs, we are primarily interested 
in the more numerous sdB stars. As mentioned above, sdB stars are understood to be relatively 
homogeneous and probably have a common evolution history from the zero age extended horizontal 
branch (ZAEHB) and beyond, while sdO stars likely follow multiple evolutionary pathways and 
might be expected to be less homogeneous and to have less simply explained properties. In the 
following sections, while we focus on the sdB stars, we will take the sdO stars "along for the ride" 
and comment on the two groups either separately (sdB or sdO), or combined as "hot subdwarfs" 
in general (both sdB and sdO). We further discuss the adopted classifications of objects from KHD 
and our treatment of them in 52.4. 



2.1. Coordinates 

To make a comparison of the KHD data with 2MASS, accurate coordinates on a consistent 
system are required. For each entry we verified the object's position by referring whenever possible 
to original published finding charts or by contacting knowledgeable observers, then locating the 
object on a chart prepared from the USNO A2.0 Catalog 4 (Monet et al. 1996). We adopted 
coordinates from USNO A2.0, except for bright objects missing from that source, in which case 



4 http://www. nofs.navy.mil/data/fchpix/ 
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we used coordinates from the Tycho-2 (H0g et al. 2000) or Hipparcos (Perryman et al. 1997; ESA 
1997) catalogs. These updated ICRS 2000.0 coordinates are accurate to a few arcseconds (~ 2" or 
better) at the epoch of observation; in many cases this represents a very substantial improvement 
over the coordinates reported by KHD. During the time in which we collected the 2MASS data, we 
had confirmed and improved coordinates for 1459 objects (not counting 11 duplicate entries), while 
for 57 objects we had not yet positively identified the correct star. These remaining 57 unrecovered 
objects are not included in the present study. The full collection of improved coordinates will be 
published elsewhere (Wade & Stark, in preparation). 

2.2. Infrared Photometry 

Infrared photometry (J, H, and K$) from the 2MASS 2IDR (which covers ~45% of the sky) 
was collected by performing searches within ~ 10" of the the updated coordinates. Visual verifica- 
tion was performed for crowded fields by comparing images from the Digitized Sky Survey 5 (DSS) 
with images from 2MASS (most hot subdwarfs are uncrowded, with no other comparably bright 
field stars within ~ 30"). Objects are included in 2MASS if they have photometric uncertainty 
< 0.155 mag (signal-to-noise ratio, S/N>7) in at least one of the three bands (J, H, or Kg)- This 
results in reliable photometry for point sources brighter than J = 15.8, H = 15.1, and K$ = 14. 3 6 . 
J — Kg errors for our objects are typically less than ~0.2 mag, with an average J — K$ error of 
~0.12 mag for our sample. Of the initial sample (1459 objects), we were able to retrieve 2MASS 
photometry for 571 hot subdwarfs (593 before removing 22 objects classified as non-subdwarfs). Of 
the 571 hot subdwarfs, only 385 had both J and K$ magnitudes listed. Of the other 186 objects, 
180 have measurements of J but only 95% confidence upper limits on K$, and six have detections 
of a nearby companion to the subdwarf (but not the subdwarf itself). Also in the observations 
returned from 2MASS, a small fraction (few percent) of subdwarfs have two possible identifications 
within ~ 10 — 20" of the input coordinates. The hot subdwarf members of these visual doubles 
were included in the present study since both finding charts and 2MASS images were available to 
identify the correct member of the pair. These visual doubles and their interpretation as possible 
physical binaries will be discussed in more detail in future work. 

Table 1 lists the identifiers and coordinates for those objects from KHD that were found in 
2MASS. Columns 1 and 2 list a sequence number 7 and the first object identifier given in KHD. 
Columns 3 and 4 contain the J2000 coordinates based on the USNO A2.0 system. Column 5 is 



5 http: //archive. stsci.edu/dss/index. html 

6 However, according to the Explanatory Supplement to the 2MASS Second Incremental Data Release (located at: 
http://www.ipac.caltech.edu/2mass/releases/second/doc/explsup.html), at high Galactic latitude, the 2MASS Point 
Source Catalog contains accurate detections 0.5 — 1.0 magnitudes fainter than these limits. 

7 The sequence number, assigned by us, is not intended to be used as a new designator in the literature, but is 
useful in navigating the KHD catalog and our Tables 1 and 2. 
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the adopted classification for the hot subdwarf (see §2.4). Column 6 contains notes on individual 
objects as well as cross references for duplicate entries (see Table 1 footnotes). Column 7 lists an 
additional identifier for the hot subdwarf, if needed, using nomenclature consistent with SIMBAD 8 . 
For example, the designator CSHK 22876—35 used by KHD for star #8 in this table is equivalent to 
the SIMBAD designator BPS CS 22876-35, and the designator KPD 00330+5229 used by KHD for 
star #55 in this table is equivalent to the SIMBAD designator KPD 0033+5229. PG 23595+1942 
used by KHD for star #1527 in this table is equivalent to the SIMBAD designator PG 2359+197. 
From now on, when we refer to a specific object from the sample, we will use both its sequence 
number (for ease in cross-referencing our Tables 1 and 2 with each other and with KHD) and a 
SIMBAD-compliant designator. If no other ID is given in the last column, the identifier from KHD 
(column 2) is consistent with SIMBAD. In this way, one SIMBAD-compliant designator for each 
star is available for cross-reference. The exceptions are the "HS" star #390 HS 1000+4704 (Heber, 
Jordan, & Weidemann 1991), the "SC" star #1165 SC 1721-3338 from Reid et al. (1988), and the 
"LBQS" stars from Berg et al. (1992, #41, 565, and 1405) which are not presently recognized by 
SIMBAD; we note that SIMBAD does use the LBQS designator for quasars from Hewett, Foltz, & 
Chaffee (1995). 



2.3. Visual Photometry 

Visual photometry (BV or uvby) for the hot subdwarfs was taken from KHD, Allard et al. 
(1994), Landolt (1992), Mooney et al. (2000), Williams, McGraw, & Grashuis (2001), and/or the 
Tycho-1 Catalog (H0g et al. 1998). For the photometry taken from KHD, when more than one set of 
photoelectric measurements was given, we favored Johnson BVover Stromgren uvby, unless the only 
B V measurement was reported as "uncertain" (then uvby was used). If only a single "uncertain" BV 
measurement was available for an object (and no Stromgren data were found), then that uncertain 
measurement was nevertheless used and treated as exact. For consistency, we took the first of these 
measurements listed in KHD, unless it was uncertain, and even if additional measurements differed 
significantly from the first. Those hot subdwarfs with only photographic visual photometry were 
not included in the optical- infrared color space examination of hot subdwarfs discussed in §3.1-3.2 
(however, those with 2MASS 2IDR detections are included in the objects discussed in §3.3). For 
stars listed in Allard et al. (1994), Landolt (1992), or Williams et al. (2001) the magnitudes from 
those sources were used in place of any reported in KHD. Stromgren uvby was taken from Mooney 
et al. (2000) for two objects: #320 PG 09146+0004 and #1327 PG 21325+1235 in preference to 
uncertain V and no BV (respectively) in KHD. Tycho-1 BV magnitudes (H0g et al. 1998; ESA 
1997) were taken for #234 BD +34° 1543, which had only a single V magnitude listed in KHD. 
No further literature search was undertaken to locate additional visual photometry. No errors are 
reported for the photometry given in KHD and we did not refer to original references to locate 



http://vizier.u-strasbg.fr/, operated at CDS, Strasbourg, France 
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errors for each individual observation. However typical errors are < 0.02 mag for those surveys 
comprising the majority of the visual photometry (i.e., Allard et al. 1994; Green et al. 1986; Mooney 
et al. 2000; Wesemael et al. 1992; Williams et al. 2001) — thus, for simplicity, all measurements of 
visual photometry were treated as exact. 

When only Stromgren photometry was available (or was taken in place of uncertain BV), we 
assumed that y = V, and b—y was converted to B — V using: 



where: mi = (v—b) — (b—y) (Turner 1990). Two objects have only synthetic b and y reported: #1076 
Balloon 083500007 and #1095 Balloon 083600004. For these two objects, B — V was estimated by 
setting mi = in this equation. 

Tycho-1 B and V magnitudes for one object, #234 BD +34° 1543, were converted to Johnson 
B and V magnitudes using: 



This conversion is valid over — 0.2 < (Bt — Vt) < 1-8, and generally gives errors < 0.015 in Vj and 
< 0.05 in (B-V)j (ESA 1997). 

Of the 571 reported hot subdwarfs found in the 2MASS 2IDR, we were able to collect pho- 
toelectric For y magnitudes for 348 stars, and B — V ox the b — y equivalent for 339 of these 348 
stars. 

All collected IR and visual photometry is presented in Table 2. Columns 1 and 2 list the 
sequence number and first identifier given in KHD. Columns 3-8 contain the 2MASS magnitudes 
(J, H, and K$) and their respective la errors. If no error is listed, then the star was not detected 
in that band, and the quoted magnitude is a 95% confidence upper limit. Columns 9-10 are the 
calculated J—K$ and combined la error for it. If no error is listed then the reported J—Ks value is 
the maximum color index when an upper limit on K$ is reported, and it is the minimum color index 
for an upper limit on J. Column 11 lists the V, y, or photographic visual magnitude of the hot 
subdwarf. Column 12 contains the B — V color index for the hot subdwarf (derived from Stromgren 
uvby or Tycho-1 BV in necessary, as discussed above). Column 13 is the adopted classification 
for the hot subdwarf (see §2.4). Column 14 contains additional notes, which are explained in the 
footnotes to Table 2. 



(B-V) = 1.584 x (b-y) + 0.681 x m x - 0.116 



Vj 

(B-V)j 



V T - 0.090 x (B T -V T ) 
0.850 x (Bt-Vt) 



2.4. 



Classifications 



During our analysis, we divided the total sample of hot subdwarfs with 2MASS 2IDR data into 
two sub-samples: "sdB/sd" and "sdO/sdOB" based on their classifications in KHD. For consistency, 
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we took the first classification given, even if there were additional conflicting classifications. Our 
"sdB/sd" group comprises those objects classified as sdB, sdB-O, and sd. Our "sdO/sdOB" group 
comprises sdO, sdO(A), sdO(B), sdO(C), sdO(D), and sdOB. There is ongoing discussion in the 
community regarding appropriate multi-dimensional classification schemes for hot subdwarf stars, 
which likely will depend on the spectral resolution used (Jeffery et al. 1997; Drilling et al. 2000, 
2003). Meanwhile, we have made no attempt to improve upon, regularize, or rationalize the various 
classifications summarized in KHD. 

As noted earlier, 22 additional objects were excluded from our sample because they are not 
subdwarfs. These objects are classified: HBB (three objects; Tables 1 and 2 sequence #: 350, 642, 
and 662), DA (two; #265 and 588), non-sd (nine; #269, 287, 574, 676, 734, 794, 1138, 1443, and 
1501), only as "bin?" (three; #732, 1076, and 1095), or had no classification (four; #398, 702, 884, 
and 967) in KHD. Lastly, #276 KUV 08388+4029 is listed as "sdB" in KHD, however, SIMBAD 
lists this object as a Seyfert 1 galaxy at z ~ 0.15. We verified its AGN nature with a low-resolution 
spectrum obtained from the Hobby-Eberly Telescope 9 , which clearly showed the Ha and H/3 broad 
emission lines at wavelengths consistent with a redshift z ~ 0.15. 

Hereafter, when we refer to "sdB stars" we mean our "sdB/sd" group, when we refer to "sdO 
stars" we mean our "sdO/sdOB" group, and when we refer to "hot subdwarfs" we mean the entire 
sample (sdB and sdO) excluding the non-subdwarfs just mentioned. A summary of the numbers 
and types of subdwarfs at each stage in the data gathering process is presented in Table 3. 

3. Hot Subdwarfs in Color Space 

3.1. Optical-IR Color-Color Plots 

After examining various combinations of color indices, we chose to focus on (J—Ks) vs. (B — V) 
and (J — Ks) vs. (V — Ks). These two combinations allow us to extract most of the information 
from the available photometry. Substituting J—H in place of J—Ks shows the same features and 
trends. Likewise, substituting V— J for V —K$ shows the same features and trends. Using H—Ks 
did not provide any additional information. 

The locations of single and composite hot subdwarfs in color space were determined by selecting 
samples of hot subdwarfs reported as single and composite from the literature. A total of 69 
reportedly single stars and 37 reported composites that had 2MASS detections were taken from the 
literature (KHD; Ferguson et al. 1984; Beers et al. 1992; Allard et al. 1994; Thejll et al. 1995; Ulla 
& Thejll 1998) or a list of hot subdwarfs spectroscopically observed and kindly provided to us by 



9 The Hobby-Eberly Telescope is operated by McDonald Observatory on behalf of The University of Texas at 
Austin, the Pennsylvania State University, Stanford University, Ludwig-Maximilians-Universitat Munchen, and 
Georg-August-Universitat Gottingen. The Marcario Low Resolution Spectrograph is a joint project of the Hobby- 
Eberly Telescope partnership and the Instituto de Astronomfa de la Universidad Nacional Autonoma de Mexico. 
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E. M. Green. After removing from this list stars with incomplete photometry and stars classified 
as "non-sd" by KHD, 33 (24 sdB, 9 sdO) of the reported composites and 59 (48 sdB, 11 sdO) of 
the reported singles remained and are used here to determine the composite and single loci (these 
stars are flagged in Table 2 with a "s" or "c" in the last column). 

When plotted in (J — Ks, B — V) or (J — Ks, V — Ks) color spaces relative to the location of 
the Pop I main sequence (Figure 1), the reported single hot subdwarfs (both sdB and sdO) occupy 
the same area in color space as other O and B stars, with the sdOs in general having slightly 
bluer colors than the sdBs. Most of the single hot subdwarfs occupy the region: (B — V) < +0.1, 
(V — Kg) < +0.2, and (J — Ks) < +0.05. This region is marked by a box on Figure 1. On the 
other hand, reported composite hot subdwarfs (both sdB and sdO) occupy a separate and distinct 
region of color space; they have redder colors than single subdwarf stars. Composite subdwarfs 
show a consistent and quite large color difference from the "typical" single subdwarf in J — Ks- 
The difference in V — Ks is also large in the mean, but the distributions of single and composite 
systems overlap slightly. The least average difference is in B — V (Figure 1). 

These patterns can be understood if the hot subdwarfs companion star in the composite 
systems is cool and relatively faint at visible wavelengths. The bluest color index, B — V, will 
be relatively unaffected by the additional light from the companion. The visual-infrared index, 
V— Ks, will betray a large infrared "excess" compared to a single hot subdwarf, but there will be a 
large spread because the index depends on the competing light from the two stars at V. The pure 
infrared index, J—Ks, will be dominated by the cool companion, whose color is very different from 
that of the hot star. 

Three reported composite sdOs have distinctly different colors in Figure 1 compared with all 
of the other reported composites. The first (labeled "1" in Figure 1) is #1355 BD — 3°5357 (or 
FF Aqr, with B — V = +0.895, J-K s = +0.752, and V-K s = +2.838), an eclipsing binary (orbital 
period ~9.2 days) classified sdOB+G8III (Dworetsky et al. 1977; Etzel et al. 1977). The second 
("2") is #1457 BD -7°5977 (with B — V = +0.541, J-K s = +0.569, and V-K s = +2.129), 
classified sdOB+K0IV-III by Viton et al. (1991). The third ("3") is #127 PG 0205+134 (with 
B-V = +0.006, J-K s = +0.847, and V-K s = +2.771), classified "sdOB comp" by Green et al. 
(1986), and further refined to sdOB+M3.5(IV?) by Allard et al. (1994) and Williams et al. (2001). 
All of these cool companions are atypically luminous for sdO-type composites, which otherwise have 
main sequence or subgiant companions reported, so it is not unreasonable for these three objects 
to have exceptionally red colors. 

The two sdB stars labeled "4" (#822 PG 1439-013) and "5" (#744 TON 183) in Figure 1, 
have anomalously blue J — Ks colors, J — Ks = —0.496 ± 0.148 and —0.494 ± 0.090 respectively 
(la errors). These colors could be due to photometric errors, or they could be intrinsic to the 
objects themselves. The color index J—H generally shows the same trend as J—Ks, and these two 
stars have essentially the same J—H colors as all other single sdBs (J—H = —0.097 and —0.196). 
Possibly, these unusual J—Ks colors are due to photometric errors in Ks- 
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One reported single sdB star, #483 TON 64 (labeled "6" in the right plot of Figure 1), has a 
rather red value for V — Kg (+0.592), yet appears "normal" in the other colors (B — V = —0.220, 
J — Kg = —0.028). Its unusual position may be due to an error in the standardization of the V 
magnitude. The two measurements for the visual magnitude of this star in KHD are discrepant: 
V = 15.86 and y = 14.59. Our adopted procedure is to favor V over y; if however we were to use 
y = 14.59 in place of V = 15.86, we would have V—Kg = —0.678 (again assuming y—Kg = V— Kg), 
which places this star near the center of the single subdwarf locus. Nevertheless for consistency, 
we continue to use the value of V for this object. 

Now placing all of the 277 hot subdwarfs with V (273 with B — V) and J and Kg measurements 
in color space (Figure 2), we observe that they occupy the same areas in color space as defined 
by the smaller sample of known composites and singles. (Objects with only upper limits on Kg 
from 2MASS are not plotted in Figure 2.) There is some additional scatter, which is most likely 
due to errors in the photometry. This scatter could, in principle, also be due in part to differences 
in interstellar extinction and reddening from object to object. Additionally, some hot subdwarfs 
are known to pulsate such as the V361 Hya stars 10 (O'Donoghue et al. 1999) and the "PG 1716" 
stars (Green et al. 2003). In the case of the V361 Hya stars, the pulsation amplitudes are generally 
less than 0.01 mag and typically only a few milli-mags with typical periods of 100-250 seconds 
(i.e., Kilkenny et al. 1997; Ostensen et al. 2001; Green et al. 2003). The "PG 1716" stars have 
peak-to-peak amplitudes < 0.05 magnitudes over timescales of about an hour (Green et al. 2003). 
Light variations due to both of these pulsations are typically less than other sources of error. 

Figure 2 suggests that candidates for composite systems (hot subdwarf+late type star) can 
be identified and targeted for future study using B — V ox V — Kg in combination with J — Kg. 
Unfortunately, in order to use this technique, accurate visual photometry is needed in addition to 
IR photometry. 

3.2. Modeling the Distribution in Optical-IR Color-Color Space 

To test the interpretation offered above, that the colors of single and composite hot subdwarfs 
may be understood if the companions are cool stars, we attempted to reproduce the location of 
composite subdwarfs in color space by combining the colors of a typical sdB with those of late type 
Pop I main sequence stars. Color indices for late type stars were taken from Johnson (1966); for 
this exercise we assumed K ss Kg. The "typical" colors for sdBs were taken from the center of 
the clump of reported single stars in Figure 1. We used {B — V) s ab ~ —0.25, (V— Kg) sdB ~ —0.75, 
(J — -Ks)sdB ~ —0.15. Based on the spread seen in Figure 1, however, it is likely that single sdB 
stars in fact have spreads in color index on the order of A(B — V) ~ ±0.1, A(V— Kg) ~ ±0.2, and 
A(J — Kg) ~ ±0.15, as might be expected from variations in the physical properties of the sdBs 



also known unofficially as "EC 14026" stars 



- 11 - 



(most notably the effective temperature); some of the spread in V— K$ and J— Kg is however due 
to limited IR photometric precision. 

The colors were combined, varying the fraction of the light at V that arises from the companion. 
The resultant locus passes among the known composite sdBs with contributions from the companion 
of ~10-60% at V, and spectral types G-K (Figure 3). These fractional contributions and companion 
types agree with those reported by Allard et al. (1994) for sdBs, and Ferguson et al. (1984) and 
Orosz, Wade, & Harlow (1997) for hot sdBs and sdOs. While this general analysis also fits the 
sdOs, the unreddened colors of a "typical" sdO should be somewhat bluer than those of a sdB, 
since sdOs are hotter than sdBs. 

The loci are plotted with the entire subdwarf sample in Figure 4. 

By introducing assumed absolute magnitudes along with colors, we can synthesize composite 
colors specifically for sdB+late type main sequence stars. We took colors and My for Pop I main 
sequence stars from Johnson (1966) and Drilling & Landolt (2000) and assumed colors and My S( jB 
for the sdB component. The sdB colors used in the previous calculation were again adopted. The 
absolute magnitudes of sdB stars are poorly constrained, with reported values spanning the range 
from My sdB ~ 3.5-6.2 (Allard et al. 1994; Bixler et al. 1991; Ferguson et al. 1984; Liebert, Saffer, 
& Green 1994; Moehler, Heber, & de Boer 1990; Orosz et al. 1997; Thejll et al. 1997). We therefore 
treated My s dB as a parameter, and varied it (choosing My s dB = 3.5, 4.0, 4.5, 5.0, and 5.5) to see 
its effect on the location of the synthesized composites in color space (Figure 5). We found that the 
locus of synthetic composites passes among the observed composite systems and best characterizes 
the observed companions as having Pop I spectral types G-K when My s dB ~ 4.5 — 5.0 mag (if 
the companions lie on the main sequence). It should be noted, however, that the same pattern 
of color loci can be created by fixing My s dB and instead varying the absolute magnitudes of the 
companions from their main sequence values. If G-K type companions were subgiants, they would 
be more luminous by perhaps ~2 magnitudes, and would therefore give the same pattern with 
a more luminous sdB, My s dB ~ 2.5 — 3.0 mag (but we would not expect to see evolved M-type 
companions). In nature, there is likely to be some range in absolute magnitude for both the sdB 
stars and their companions. (In addition to the variations of luminosity and bolometric correction 
along the zero-age EHB, sdB stars should evolve to higher luminosity on timescales of 10 8 years.) 

The (main-sequence) loci are plotted with the entire hot subdwarf sample in Figure 6. 

These simple models give results consistent with previous studies concerning the companion 
spectral types and their light contributions. We cannot, however, assume an My for the companion 
and invert the calculation case-by-case to give reliable individual values for My S( jb, nor vice versa, 
because of the uncertainties in the measured magnitudes, colors, and reddening. Additionally there 
is the expected intrinsic spread in My s dB referred to above (for further discussion see Saffer et al. 
1994). For similar reasons, neither can we infer exact companion spectral types for individual 
objects. 

Of the 33 currently known composite hot subdwarf systems, essentially all of the cool com- 
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panion stars for which a spectral class is known are of spectral types G or K. Our models show 
that composite systems with G, K, and early M-type companions show the most significant differ- 
ences in color from a single sdB (particularly in J—Ks) and from main sequence stars. This could 
explain why most known composite systems have GK-type companions: for companions earlier 
then ~G, the combined system looks more like a main sequence star (or a cool, Pop n subdwarf) 
than a sdB star. Such systems would likely be excluded from catalogs of hot subdwarfs. For 
main-sequence companions much later than ~M0, the effect on the combined system light is so 
small that the companion would not easily be detected in the visual. Most of the spectroscopic 
studies of composite sdBs have been done in the visual and blue region of the spectrum where 
very little color difference between a composite sdB+dM and a single sdB exists, and there are 
indeed few identifications of dM companions — thus far there are only four known: PG 1241—084 
(also HW Vir; Menzies & Marang 1986), PG 1336-018 (Kilkenny et al. 1998), HS 0705+6700 
(Drechsel et al. 2001), and PG 1017-086 (Maxted et al. 2002). All of these are very short period 
binaries (2 h 48 m , 2 h 25 m , 2 h 18 m , and l h 45 m respectively), and the the nature of the companion was 
studied in part by eclipses or reflection effects in the light curve. PG 1241—084 (KHD sequence 
#605) and PG 1017-086 (#407) are found in the 2MASS 2IDR. Colors for such close binary stars, 
constructed from non-simultaneous measurements of visible and infrared magnitudes, may be un- 
usual. The color indices for these two stars are nevertheless consistent with single hot subdwarfs 
(#605 PG 1241-084: J-K s = -0.149, B-V = -0.332, V-K s = -0.110; #407 PG 1017-086: 
J-K s = -0.176, B-V = -0.225, V-K s = -0.417). A few other objects can be noted in the 
region of color space that, based on our simple models, would be occupied by sd+dM composites 
(Figures 4 and 6). Owing to uncertainties in the measured magnitudes and colors, however, we can- 
not clearly distinguish early M companions from late K companions — near-IR or IR spectroscopy 
will be needed to distinguish clearly between late K and early M-type companions. Yet it seems 
clear, based on our color plots, that there are not great numbers of "hidden" early M-type (~M0) 
companions to hot subdwarfs. For main sequence companions substantially later than ~M0, the 
J—Ks color is rapidly returning to that of a single subdwarf. 

Identifying and characterizing further M-type companions to sdB stars would help settle some 
of the questions about the nature of earlier-type (G-K) companions, in particular whether they are 
predominately main sequence or subgiant stars. Late M-type companions can safely be assumed to 
be on the main sequence. If many of these are found among the sdB stars with "composite" colors, 
the relatively faint My^dB inferred from such systems would imply that the G-K type companions 
in other composite-color systems are also close to the main sequence. If late-M companions are 
not found in great numbers, it could be inferred that the G-K companions are mostly subgiants, 
since the corresponding bright My jS dB inferred from this assumption would explain the absence of 
sd+dM systems detectable from photometry. The existence of many subgiant companions, with 
their short lifetimes as subgiants, would be interesting, since the evolutionary timescale of sdB stars 
is also short. 
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3.3. Hot Subdwarf Distribution in J—Ks Color 

Since visual photometry is not available for all of the hot subdwarfs contained in the 2IDR, and 
since there is a more distinct separation between composites and single subdwarfs stars in J — Ks 
than in the other colors, we examined the distribution of subdwarfs based on their J — Ks colors 
alone. A histogram of J—Ks for all subdwarfs with J and K$ measurements (385 subdwarfs: 294 
sdB, 91 sdO), reveals a bimodal distribution (Figure 7). One peak is centered near J—Ks ~ —0.15 
(~B2-B3 spectral type), and the other is centered near J—Ks ~ +0.30 (~G5 spectral type). There 
are also wings stretching to both redder and bluer colors. This distribution cannot be explained 
by reddening, because no significant trends were observed with galactic latitude. We would also 
expect reddening to produce a tail to the red from a single peak, but not a pronounced second 
peak. Finally, the amount of extinction that would be required is also excessive. 

We sorted the sdB and sdO samples into three groups by J magnitude. The bimodal peaks 
are very pronounced in the brightest third of each sample, less distinct in the middle third, and 
practically non-discernible in the faintest third. This pattern can be fully explained by the increase 
in photometric uncertainty at fainter magnitudes. 

In keeping with the results found from our modeling of color-color distributions, we postulate 
that the stars within the blue peak of the J—Ks distribution {J—Ks < +0.05) are single stars (or 
sd+WD pairs, with colors indistinguishable from single hot subdwarfs), and that stars in the red 
peak (J — Ks > +0.05) are composite (sd+late type) systems. The color index separating these 
two groups of stars was adopted from the box defining the region of single subdwarfs in Figure 1. 
With this definition, we find that 58+5% of the subdwarf stars in our 2MASS 2IDR sample are 
apparently single, and 42+5% are composite. The quoted 95% confidence errors are sampling errors 
computed from the binomial distribution, taking the size of the sample into account. If the color 
index of the boundary separating single and composites is increased (decreased) by 0.05 mag, the 
fraction of subdwarfs that are "blue" (i.e., single) is increased (decreased) by 3% (3%). 

While the 2IDR sample is representative of the KHD catalog as a whole, it is likely that neither 
of these is representative of the true hot subdwarf population in the solar neighborhood. There are 
the usual completeness issues, including heterogeneous magnitude limits and different classification 
criteria used by the sources from which KHD was compiled (although we note that about 2/3 of the 
cataloged subdwarfs are drawn from the PG survey, Green et al. 1986, which has rather well-defined 
properties). We draw renewed attention, however, to a kind of censorship bias that is present in the 
PG survey. An unknown number of composite subdwarfs may have been excluded "by design" (i.e., 
not included in the PG catalog and never identified as hot subdwarfs), since according to Green 
et al. (1986) UV excess stars were omitted from their catalog if their spectra showed a pronounced 
Ca II K line. (The supposition was that many of these would be metal-weak F or G subdwarfs, which 
met the purely photometric criterion for inclusion as UV excess stars owing to a combination of 
reduced line blanketing and larger than average photometric errors.) Thus composite hot subdwarf 
systems with F or G companions might be underrepresented in KHD and the present sub-sample 
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drawn from it. (On the other hand, since some PG subdwarfs are found in other independent 
surveys, some of the composite systems "missing" from PG may nevertheless have been included in 
the KHD compilation, from other sources.) The number of "excluded" objects is roughly 1000-1200 
(R. F. Green 2003, private communication), but the fraction of these that may be composite hot 
subdwarfs rather than G subdwarfs is not known, so these selection effects are difficult to quantify. 
In time, objective prism-based surveys such as the one being carried out with the Hamburg Schmidt 
telescope may provide a catalog of hot subdwarfs that is more clearly representative of the Galactic 
population. Meanwhile, we proceed to analyze the 2MASS 2IDR sample from the KHD catalog. 
It must be realized, however, that since the 2MASS 2IDR represents a magnitude limited survey, 
the composite systems (being more luminous than single subdwarfs) will be overrepresented in this 
sample, other things being equal (see §4 for further discussion of this point). 

The same bimodal distribution in J — K$ exists in the sub-sample of sdB/sd stars that is 
found in the full sample (58±6% have J-K s < +0.05, and 42±6% have J-K s > +0.05). The 
sub-sample of sdO/sdOB stars shows a less clear pattern. There is a hint of three groups: one near 
J — Ks ~ —0.2, another near J — Ks ~ +0.2, and the third near J — Ks ~ +0.5, but the small 
number of stars in each group makes this finding quite tentative. The two larger groups of sdO 
stars (J—Ks ~ —0.2 and +0.2), if real, can be considered to fall to the "blue" side (lower J—Ks) of 
the the two peaks in the sdB distribution, consistent with the higher temperatures of the sdO with 
respect to the sdB stars. Overall, the observed bimodal distribution in J—Ks of all hot subdwarfs 
is defined by the majority sdB stars. The fraction in each peak of the J—Ks distribution does not 
significantly change if the sdO stars are removed. A breakdown of the fractions of subdwarfs by 
J — Ks color is shown in Table 4. 



3.4. Fitting the sdB J—Ks Distribution 

We have fitted the bimodal J — Ks distribution of the sdB stars (alone) to the sum of two 
Gaussians, centered on the blue and red peaks. An initial "fit-by-eye" to the binned data was 
improved by minimizing x 2 on a range that was trimmed (—0.45 < J—Ks < +0.6) to remove the 
effects of the reddest and bluest outliers. The x 2 minimization was done using a systematic gridded 
search in the neighborhood of the initial fit. There are six fitting parameters: two mean values 
(or centers), two dispersions (or widths), and two normalizations (or amplitudes). The final fit has 
X 2 = 6.8 (reduced Xr = 0.45, with 15 degrees of freedom). The parameters for the fit are listed in 
Table 5 and the fit is shown in Figure 8. 

While the x 2 statistic provides a useful criterion for improving the fit, its final numerical value 
should be used with caution, since the range of the fit was trimmed, and since some of the data 
bins contain only a few members. In particular, the 95% confidence intervals shown in Table 5 are 
indicative rather than definitive. The small value of X\\ nevertheless demonstrates the plausibility 
of a two-Gaussian fit, which is confirmed by visual inspection. The "blue" Gaussian has a somewhat 
smaller Gaussian dispersion than the red one, but both dispersions (a « 0.09, 0.16) are close to the 
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average photometric error for this sample from 2MASS, a{J — Ks) ~ 0.12. This makes it difficult 
to state whether there is any intrinsic spread in the J—Ks colors of the single or composite sdBs. 

To the extent that the separation of single and composite sdB stars in the sample is accurately 
represented by two Gaussians, it is seen that there is very little overlap in J—Ks color index between 
the two groups of stars. The areas of the two Gaussians are in the proportion 57:43, consistent 
with the binomial result found using a sharp demarcation at J —Kg = 0.05. 

4. A "Volume Limited" Hot Subdwarf Sample 

4.1. Defining the Sample 

There is a strong advantage in working with the J — Ks color index to discuss the fraction 
of composite-color binaries among hot subdwarfs, as in §3.3 and 3.4 above, since the data are 
available for a larger sample than would be the case for optical or mixed optical-infrared colors. 
Also, the J and K$ from 2MASS are homogeneous in nature. However, a strong selection bias in 
these data exists, because 2MASS itself is an approximately magnitude limited survey (albeit with 
different magnitude limits in each of three passbands). The entire sample of hot subdwarfs we have 
been discussing (common to both KHD and 2MASS 2IDR) is thus a magnitude limited sample 
in some sense. The volume sampled for single subdwarfs is therefore smaller than that sampled 
for composites, since the presence of the late type companion increases the distance at which a 
composite system remains brighter than the limiting 2MASS magnitude. 

To further explore this selection bias, we constructed (J, J — Ks) and (Kg, J — Ks) color- 
magnitude diagrams for the full 2MASS 2IDR sample (Figure 9). Several interesting patterns are 
noted. First, the two populations (single and composite) can be seen, and both are observed to 
cover the entire magnitude range (i.e., composites and/or singles are not "clumped" at a specific 
apparent magnitude). Second, a majority of the objects with extreme J—Ks colors (either blue or 
red) are found at the faintest magnitudes, an effect largely explained by photometric errors. Third, 
it can be seen that the sample is limited by detections in the K$ band. There is a clear cut-off at 
Ks ~ 15.75, while in J the cut-off magnitude varies by color (redder objects that are detected in 
both J and K$ run to a fainter J magnitude than the bluer objects). 

The interpretation of this third feature is that a blue (single) object detected near J = 15.5 
or fainter is missing at Kg and thus is not part of the plotted sample; a red (composite) object at 
the same J will also have a K$ detection and thus will appear in Figure 9. (This effect is distinct 
from the fact that composites appear brighter at J due to the combined light from two stars.) 

Recognizing this selection bias, we can attempt to define an approximately volume limited 
sample of hot subdwarfs. (To our knowledge, this is the first time that this bias has been discussed in 
the context of hot subdwarf composites. Theissen et al. (1992) discussed a "statistically complete" 
sample containing 11 hot subdwarfs, however they did not address the composite vs. single bias.) 
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To define our volume limited sample, we try to remove those subdwarfs that we observe to be 
composite, but which would not have been detected by 2MASS had they been single. With only 
a single color index to work with, and with uncertainties in both the photometry and the relative 
luminosities of the hot subdwarfs and their companions, we cannot rely on such a technique to say 
whether any given individual system lies within the defined volume. The sample so constructed 
can be used in a statistical sense, however, to reflect the general nature of the true distribution, 
and at least to indicate the size of the selection effect. 

The observed cut-off magnitude from Figure 9 for a single hot subdwarf is Ks ~ 15.75. Taking 
this limit as the apparent magnitude of the hot subdwarf component in a composite system, we 
can compute the corresponding apparent magnitude for the total light from the system, using the 
models discussed in §3.2 and shown in Figures 3 and 4. This total magnitude will vary as the 
spectral type and/or relative luminosity of the late type companion is varied. The locus of thus- 
computed magnitudes is shown in the lower panel of Figure 9 for a hot subdwarf with the "typical" 
color J — Ks = —0.15 adopted in §3.2, for companions with colors of the Pop I main sequence 
spectral types GO, K0, K5, and MO. Along each of these lines, the color and total magnitude varies 
as the relative brightnesses of the two component stars are varied. Note that no assumption has 
been made about the absolute magnitude of either the hot subdwarf or the cool companion. 

We now replace the family of lines (corresponding to different spectral types of the cool star) 
with a single representative line. (This is because based on J — Ks alone we cannot know the 
true spectral types of the late type companions without making assumptions about the relative 
luminosities of the component stars.) We use this line to approximate the appropriate magnitude 
limit for a volume limited sample, again in a statistical sense. The line is chosen to pass though the 
middle of the known composite sdBs shown in Figure 3; it thus represents a ~ 5% contribution (at 
V) by a K5 star, or a contribution of 40-80% at K0 in Figure 9. The change in the volume-limiting 
Kg magnitude is thus most rapid near J — Ks ~ +0.3 (see Figure 3). All hot subdwarf stars in 
the (Ks, J — Ks) color-magnitude diagram that fall above this line are included in our statistical 
volume limited sample. The resulting sample size is 300 hot subdwarfs (228 sdB and 72 sdO). The 
average error in J — Ks for this new sample is 0.10 mag (slightly reduced from 0.12 mag, which 
applies to the full J—Ks sample). 



4.2. Color-Color and J —Ks Distribution of the Volume Limited Sample 

For the volume limited sample, visual magnitudes (B For the uvby equivalent) are available 
for 176 sdB stars (4 additional with V only) and 61 sdO stars (1 additional with V only). The 
volume limited sample occupies the same regions of (J — Ks, B — V) and (J — Ks, V — Ks) color 
spaces (Figure 10) that the magnitude limited sample does (compare Figures 2 and 10), but many 
of the more extreme red objects (in J —Ks) have been removed. 

The J—Ks histogram made using only the volume limited sample (Figure 11) can be compared 
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with the corresponding diagram for the entire magnitude limited sample (Figure 7). The distri- 
bution is still bimodal, with one peak centered near J — Ks ~ —0.15 as before (~B2-B3 spectral 
type), and the other near J — Ks ~ +0.30 (~G5 spectral type). The sdB stars still dominate the 
total hot subdwarf distribution, and the sdO stars show a stronger hint of bimodality (again with 
both peaks shifted slightly blueward of the sdB peaks). The blue wing to the distribution remains, 
while the red wing has been greatly reduced. Most important, the fractions of single ("blue") and 
composite ("red") systems have been altered (Table 6). Whereas the magnitude limited sample 
showed a 58:42 split between single and composite sdB stars (Table 4), the volume limited sample 
shows a 75:25 split in favor of single sdB stars. Likewise, the split for the total hot subdwarf sample 
is 73:27 in the volume limited sample (as opposed to 58:42 in the magnitude limited sample). 

We can again describe the sdB J— Kg distribution for the volume limited sample as the sum 
of two Gaussians. The new x 2 = 8.9 (reduced Xr = 0.59, with 15 degrees of freedom). The 
parameters for the fit are given in Table 7, and the fit is shown in Figure 12. 

The parameters of the blue Gaussian fits for the magnitude limited and volume limited samples 
are essentially the same. This follows because the number of stars that fall in the blue peak was 
not substantially altered in creating the volume limited sample. The width of the blue Gaussian 
is consistent with the average photometric error from 2MASS for this sample, a(J — Ks) ~ 0.10, 
implying that there is little intrinsic spread in the J—Ks color index of single sdB stars. 

The red Gaussians for the magnitude limited and volume limited samples have different mean 
values (centers) , with the volume limited sample shifted blueward relative to the magnitude limited 
sample [A(J — Kg) ~ —0.04]. The normalizations (amplitudes) are different, owing to to the 
decrease in the number of red stars in the volume limited sample. The change in center is likely 
due to the removal of the reddest and also the faintest (and thus most uncertain) measurements 
from the red peak. The dispersion of the red Gaussian that represents the volume limited sample, 
a(J—Ks) ~ 0.10, is reduced from the magnitude limited case, and is now consistent with the 2MASS 
photometric error for this new sample (« 0.10 mag); this implies that there is little variation in the 
J—Ks color of the remaining composite systems. 

5. Summary 

Using updated and improved coordinates for hot subdwarfs cataloged by KHD, we have re- 
trieved 2MASS 2IDR photometry for 571 hot subdwarfs. This is the largest sample of hot subdwarfs, 
and in particular composite systems, ever examined. We divided the sample into two sub-samples 
based on their classifications in KHD: sdB/sd ("sdB" group), sdO/sdOB ("sdO" group). We ex- 
amined these two sub-samples both independently and combined as a whole. We examined the 
distribution of hot subdwarfs in optical-IR color-color plots, by combining 2MASS with visual 
photometry (B V oi uvby) from KHD (or a few selected other sources). The most revealing combi- 
nations of color indices we found were: (J—Ks, B — V) and (J—Ks, V—Kg)- We also examined the 
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distribution of hot subdwarfs in J—Ks only and found a bimodally distributed population. Finally, 
we attempted to define an approximately volume limited sample of hot subdwarfs for statistical 
purposes. This is the first time (to our knowledge), that the selection bias favoring composite hot 
subdwarfs in magnitude limited surveys has been discussed, and an approximately volume limited 
statistical sample defined. 

In keeping with the exploratory nature of this work and the large sample size, we have favored a 
statistical approach, in contrast to the attempt to find a best characterization of individual binaries. 
(The relatively large mean errors of the 2MASS photometry also figured in the decision to use this 
approach.) In the same spirit, we have contented ourselves with "naive" modeling of the colors, e.g., 
adopting single values for the colors and magnitudes of the subdwarfs rather than a distribution of 
these quantities — our goal in the modeling is to support the overall division of the sample of stars 
into single and composite groups, by demonstrating the plausibility of such modeling. 

We found that the locations of hot subdwarfs in visual-IR color-color diagrams can be modeled 
(very simply and without assuming luminosities for either the hot subdwarf or the companion) as 
either a single B-type star or as a combination of two stars (spectral type B representing the hot 
subdwarf, and spectral type G-K for a companion) with contributions from the late type companion 
of ~ 10 — 60% at V. To a first approximation sdO stars can be described in the same manner as 
sdB stars. If main sequence companions are adopted, then the sdB stars are best described with 
Mv,sdB ~ 4.5 — 5.0 mag, and they typically have G-K type companions. However, the observed 
distribution can be reproduced equally well by assuming subgiant companions with more luminous 
sdB stars (My jS dB ~ 2.5 — 3.0 mag) — our current data cannot distinguish between these two 
possibilities. In a histogram of the IR color index J—Ks, the two peaks of the bimodal distribution 
can be understood as single stars (blue peak at J — Ks = —0.158) and composite systems (red 
peak at J — Ks = +0.311). This bimodal distribution is also present in the approximately volume 
limited sample, again with the two peaks at J — Ks = —0.154 and +0.276. Making a cut at 
J — Ks = +0.05, we found the ratio between single and composite sdB stars is 58:42 (favoring 
singles) in the magnitude limited sample, while in the volume limited sample it is 75:25 (favoring 
singles). Note that these ratios do not include any extrapolation to companions that do not affect 
the photometry (i.e., WD or late M dwarfs), since such extrapolation depends on the assumptions 
made about the absolute magnitude of the subdwarfs or their companions. The J—Ks distribution 
for both the magnitude and volume samples can be described as the sum of two Gaussians with 
very little overlap and with fractions of single vs. composite in agreement with a sharp cut at 
J—Ks = +0.05. An interpretation of these ratios is difficult due to biases in the initial subdwarf 
catalog. 

Comparison to other surveys of composite systems is difficult due to the unknown selection 
effects involved. Other surveys have generally concluded that > 50% of hot subdwarfs have late 
type companions (e.g., Allard et al. 1994; Bixler et al. 1991; Ferguson et al. 1984). However, 
these percentages depend greatly on the corrections made to account for late type companions that 
were undetectable by the survey (i.e., M-dwarfs) and the assumed luminosity of the subdwarfs or 
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observed companions (i.e., main sequence or subgiant). Based on their observations of composite 
spectrum sdB stars from the PG survey (Green et al. 1986), Ferguson et al. (1984) extrapolated 
that about 50% of the sdB stars in the PG catalog had main sequence companions of G8 or later 
(and a few had companions earlier than G8). Allard et al. (1994) identified 31 composite candidates 
out of 100 sdB stars observed photometrically, and deduced a fraction of sdB stars with cool main 
sequence companions of 54-66%, when allowance was made for perceived selection effects. From 
intermediate resolution spectroscopy of the Ca II IR triplet region of 40 sdB stars, and an assumed 
"typical" luminosity of a sdB star, Jeffery & Pollacco (1998) claim that the seven composite systems 
they observed contained subgiant companions. Therefore they conclude that the the binary fraction 
must be much higher than what Allard et al. (1994) concluded based on assuming main sequence 
companions. From a sample of 40 UV selected hot subdwarfs, Bixler et al. (1991) found that eight 
showed composite spectra, and calculated that 65-100% of hot subdwarfs had a late type main 
sequence companion. However, our volume limited sample shows significantly fewer composites 
then our magnitude limited sample. Thus, by studying only magnitude limited samples (or even 
less well-defined samples) of hot subdwarfs, it is conceivable that the total number of composite 
systems could have been significantly over-estimated by these previous surveys (all other factors 
being equal). 

The complete 2MASS catalog was released as this paper was nearing completion. We will 
address the full sky sample in a future contribution. 

The authors thank E. M. Green for providing a list of spectroscopically observed hot subd- 
warfs, G. B. Berriman for his help interpreting the 2MASS catalog, and R. F. Green for helpful 
correspondence. RW thanks STScI for its hospitality. 

This research has been supported in part by: NASA Grid Stars grant NAG5-9586, NASA 
GSRP grant NGT5-50399, and a NASA Space Grant Fellowship through the Pennsylvania Space 
Grant Consortium. 

This publication makes use of: data products from the Two Micron All Sky Survey, which is 
a joint project of the University of Massachusetts and the Infrared Processing and Analysis Cen- 
ter/California Institute of Technology, funded by the National Aeronautics and Space Administra- 
tion and the National Science Foundation; the USNOFS Image and Catalogue Archive operated by 
the United States Naval Observatory, Flagstaff Station (http://www.nofs.navy.mil/data/fchpix/); 
the SIMBAD database, operated at CDS, Strasbourg, France; and the Digitized Sky Surveys. The 
Digitized Sky Surveys were produced at the Space Telescope Science Institute under U.S. Govern- 
ment grant NAG W-2166. The images of these surveys are based on photographic data obtained 
using the Oschin Schmidt Telescope on Palomar Mountain and the UK Schmidt Telescope. The 
plates were processed into the present compressed digital form with the permission of these insti- 
tutions. 



-20- 



REFERENCES 

Allard, F., Wesemael, F., Fontaine, G., Bergeron, P., k, Lamontagne, R. 1994, AJ, 107, 1565 

Allard, F., Wesemael, F., Lamontagne, R., Bergeron, P., & Fontaine, G. 1987, in IAU, Colloquium 
no. 95, The Second Conference on Faint Blue Stars, ed. A. G. D. Philip, D. S. Hayes, & 
J. W. Liebert (Schenectady, NY: L. Davis Press Inc.), 611 

Baade, W. 1944, ApJ, 100, 137 

Beers, T. C, Preston, G. W., Shectman, S. A., Doinidis, S. P., & Griffin, K. E. 1992, AJ, 103, 267 
Berg, C, Wegner, G., Foltz, C. B., Chaffee, F. H., Jr., & Hewett, P. C. 1992, ApJS, 78, 409 
Bixler, J. V., Bowyer, S., & Laget, M. 1991, A&A, 250, 370 

Brown, T. M., Ferguson, H. C, Davidsen, A. F., & Dorman, B. 1997, ApJ, 482, 685 

Brown, T. M., Bowers, C. W., Kimble, R. A., Sweigart, A. V., & Ferguson, H. C. 2000, ApJ, 532, 
308 

Caloi, V. 1972, A&A, 20, 357 

Drilling, J. S., & Landolt, A. U. 2000, Allen's Astrophysical Quantities, 4th Edition, ed. A. Cox 
(New York: Springer- Verlag) , 381 

de Boer, K. S., Drilling, J., Jeffery, C. S., & Sion, E. M. 1997, in The Third Conference on Faint 
Blue Stars, ed. A. G. D. Philip, J. W. Liebert, & R. A. Saffer (Schenectady, NY: L. Davis 
Press Inc.), 515 

Drechsel, H., Heber, U., Napiwotzki, R., Cstensen, R., Solheim, J.-E., Johannessen, F., Schuh, 
S. L., Deetjen, J., & Zola, S. 2001, A&A, 379, 893 

Drilling, J. S., Moehler, S., Jeffery, C. S., Heber, U., & Napiwotzki, R. 2000, in The Kth Reunion, 
ed. A. G. D. Philip (Schenectady, NY: L. Davis Press), 49 

Drilling, J. S., Moehler, S., Jeffery, C. S., Heber, U., & Napiwotzki, R. 2003, in Probing the 
Personalities of Stars and Galaxies, ed. A. G. D. Philip (Schenectady, NY: L. Davis Press), 
in press 

Dworetsky, M. M., Lanning, H. H., Etzel, P. B., & Patenaude, D. J. 1977, IBVS 1284 
ESA 1997, The Hipparcos & Tycho Catalogues, ESA SP-1200 

Etzel, P. B., Lanning, H. H., Patenaude, D. J., & Dworetsky, M. M. 1977, PASP, 89, 616 
Ferguson, D. H., Green, R. F., & Liebert, J. 1984, ApJ, 287, 320 
Green, E. M., et al. 2003, ApJ, 583, 31 



- 21 - 



Green, R. F., Schmidt, M., & Liebert, J. 1986, ApJS, 61, 305 
Greenstein, J. L., & Sargent, A. I. 1974, ApJS 28, 157 
Greggio, L., & Renzini, A. 1990, ApJ, 364, 35 
Heber, U. 1986, A&A, 155, 33 

Heber, U., Jordan, S., & Weidemann, V. 1991, in White Dwarfs, Proceedings of the 7th European 
Workshop, ed. G. Vauclair & E. Sion (Norwell, MA: Kluwer Academic Publishers), 109 

Heber, U., Moehler, S., Napiwotzki, R., Thejll, P., & Green, E. M. 2002, A&A, 383, 938 

Hewett, P. C, Foltz, C. B., & Chaffee, F. H. 1995, AJ, 109, 1498 

H0g, E., et al. 1998, A&A, 335, L65 

H0g, E., et al. 2000, A&A, 355, L27 

Humason, M. L., & Zwicky, F. 1947, ApJ, 105, 85 

Jeffery, C. S., & Pollacco, D. L. 1998, MNRAS, 298, 179 

Jeffery, C. S., Drilling, J. S., Harrison, P. M., Heber, U., Moehler, S. 1997 A&AS, 125, 501 
Johnson, H. L 1966, ARA&A, 4, 193 

Kilkenny, D., Heber, U., & Drilling, J. S. 1988, SAAOC, 12, 1; and expanded electronic version 
1992, D. Kilkenny private communication (KHD) 

Kilkenny, D., Koen, C, O'Donoghue, D., & Stobie, R. S. 1997, MNRAS, 285, 640 

Kilkenny, D., O'Donoghue, D., Koen, C, & van Wyk, F. 1998, MNRAS, 296, 329 

Landolt, A. U. 1992, AJ, 104, 340 

Liebert, J., Saffer, R. A., & Green, E. M. 1994, AJ, 107, 1408 

Maxted, P. F. L., Heber, U., Marsh, T. R., & North, R. C. 2001, MNRAS, 326, 1391 

Maxted, P. F. L., Marsh, T. R., Heber, U., Morales-Rueda, L., North, R. C, & Lawson, W. A. 
2002, MNRAS, 333, 231 

Mengel, J. G., Norris, J., & Gross, P. G. 1976, ApJ, 204, 488 

Menzies, J. W., & Marang, F. 1986, IAUS, 118, 305 

Moehler, S., Heber, U., & de Boer, K. S. 1990, A&A, 239, 265 

Monet, D., et al. 1996, USNO-SA1.0, (U.S. Naval Observatory, Washington DC) 



- 22 - 



Mooney, C. J., Rolleston, W. R. J., Keenan, F. P., Pinfield, D. J., Pollacco, D. L., Dufton, P. L., 
& Katsiyannis, A. C. 2000, A&A, 357, 553 

O'Connell, R. W. 1999, ARA&A, 37, 603 

O'Donoghue, D., Koen, C, Kilkenny, D., Stobie, R. S., & Lynas-Gray, A. E. 1999, in ASP Conf. 
Ser. 169, 11th European Workshop on White Dwarfs, ed. J.-E. Solheim, & E. G. Meistas 
(San Francisco: ASP), 149 

Orosz, J. A., Wade, R. A., & Harlow, J. J. B. 1997, AJ, 114, 317 

Gstensen, R., Heber, U., Silvotti, R., Solheim, J. E., Dreizler, S., & Edelmann, H. 2001, A&A, 378, 
466 

Perryman, M. A. C, et al. 1997, A&A, 323, L49 

Reid, N., Wegner, G., Wickramasinghe, D. T., Bessell, M. S. 1988, AJ, 96, 275 
Saffer, R. A., Bergeron, P., Koester, D., & Liebert, J. 1994, ApJ, 432, 351 

Saffer, R. A., Green, E. M., & Bowers, T. P. 2000, in ASP Conf. Ser. 226, 12th European Conference 
on White Dwarfs, ed. J. L. Provencal, H. L. Shipman, J. MacDonald, & S. Goodchild (San 
Francisco: ASP), 408 

Theissen, A., Moehler, S., Heber, U., Schmidt, J. H. K., & de Boer, K. S. 1995, A&A, 298, 577 

Theissen, A., Moehler, S., & Schmidt, J. H. 1992, in Lecture Notes in Physics no. 401, The At- 
mospheres of Early- Type Stars, ed. U. Heber, & C. S. Jeffery (New York: Springer- Verlag) , 
264 

Thejll, P., Flynn, C, Williamson, R., & Saffer, R. 1997, A&A, 317, 689 
Thejll, P., Theissen, A., & Jimenez, R. 1994, A&A, 292, 457 
Thejll, P., Ulla, A., & MacDonald, J. 1995, A&A, 303, 773 
Turner, D. C. 1990, PASP, 102, 1331 
Ulla, A., & Thejll, P. 1998, A&AS, 132, 1 

Ulla, A., Zapatero Osorio, M. R., Perez Hernandez, F., & MacDonald, J. 2001, A&A, 369, 986 
Viton, M., Deleuil, M., Tobin, W., Prevot, L., & Bouchet, P. 1991, A&A, 242, 175 
Wesemael, F., Fontaine, G., Bergeron, P., Lamontagne, R., & Green, R. F. 1992, AJ, 104, 203 
Williams, T., McGraw, J. T., & Grashuis, R. 2001, PASP, 113, 490 



This preprint was prepared with the AAS IATgjX macros v5.0. 



-23- 




-0.4 0.4 0.8 -1 1 2 3 

(B-V) (V-K s ) 



Fig. 1. — Color-color diagrams for reported single (circles) and composite (triangles) hot subdwarfs. 
Filled points represent hot subdwarfs in the sdB group (48 single, 24 composite), open points 
represent stars in the sdO group (10 single, 8 composite). The left panel shows (J—Kg) vs. (B — V), 
while the right panel shows {J—K$) vs. (V—Ks). The diagonal solid line in each panel indicates the 
location of the Population I main sequence (with the approximate Pop I spectral types indicated 
along the B — V and V — Ks axes). The single hot subdwarfs are mostly contained within the box 
defined by: B — V < +0.1, V—Ks < +0.2, and J—K$ < +0.05, which is shown in each panel with 
a dotted line. Objects labeled 1-5 in both panels, and 6 in the right panel, are discussed further in 
the text (§3.1). The arrows indicate color shift due to one magnitude of extinction at V. 



-24- 




Fig. 2. — Same as Figure 1, except all hot subdwarfs with 2MASS 2IDR and visual photometry 
are shown. SdBs are denoted with circles and sdOs with crosses ("x"). There are 201 sdBs and 
71 sdOs in the left panel and 206 sdBs and 71 sdOs in the right; these totals include the reported 
composite and single stars of Figure 1. One sdB, #462 TON 56, has a published B — V = —0.90 
which places it off the left side of the left panel, this BV photometry from Green et al. (1986) is 
listed as uncertain; however, it is shown in the right panel (V — K$ = +0.263, J — K$ = +0.279). 
Four other sdBs are missing from the left panel (but not the right), because they have no B (or 
equivalent) measurements reported. 
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(B-V) (V-K s ) 

Fig. 3. — Grid of composite colors computed by combining the light from a hot subdwarf (B — V ss 
—0.25, V — Ks ~ —0.75, J — K$ ~ —0.15) with that of "typical" Pop I main sequence stars (GO, 
G5, K0, K5, and MO; Johnson 1966) assuming varying fractional contributions to the total light 
in the V band by the late type star (fractions contributed by the companion: 0, 1, 5, 10, 20, 30, 
40, 50, 60, 80, and 100%). Circles are the reported single subdwarfs (solid: sdB, open: sdO) and 
triangles are the reported composite subdwarfs (solid: sdB, open: sdO) from Figure 1. The large 
arrows indicate the color change due to one magnitude of extinction at V. 
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Fig. 4. — Same as Figure 3, except all sdBs (circles) and sdOs (x) with 2MASS 2IDR photometry 
are shown, as in Figure 2. 
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Fig. 5. — Grid of composite colors computed by combining the light from a hot subdwarf (B — V ~ 
-0.25, V-K s ~ -0.75, J-K s ~ -0.15) with that of "typical" Pop I main sequence stars (F0, F5, 
GO, G5, K0, K5, and MO; Johnson 1966) and assuming values of My for the subdwarf (My ;S( jB = 3.5, 
4.0, 4.5, 5.0, and 5.5) — the line for My, s dB = 4.5 has been darkened to help guide the eye. The 
lines for G5 and F5 have been left off the right plot for clarity. Circles are the reported single 
subdwarfs (solid: sdB, open: sdO) and triangles are the reported composite subdwarfs (solid: sdB, 
open: sdO) from Figure 1. The large arrows indicate the color change due to one magnitude of 
reddening. 
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Fig. 6. — Same as Figure 5, except all sdBs (circles) and sdOs (x) with 2MASS 2IDR photometry 
are shown as in Figure 2. 
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Fig. 7. — Histogram of J — K$ for all hot subdwarf stars found in the 2MASS 2IDR (bin size 
0.05 mag). The solid line indicates the full sample of hot subdwarfs, the dashed line indicates the 
sdB group, and the dotted line indicates the sdO group (these groups were defined in §2). Single 
subdwarfs are believed to reside in the left peak of the distribution centered around J—K$ ~ —0.15, 
while composite systems (sd+later type) reside in the right peak centered around J—K$ ~ +0.3, 
with the sdOs falling to slightly bluer colors than the sdBs in general. A tabulation of the numbers 
and fractions of single and composite hot subdwarfs can be found in Table 4. The arrow indicates 
the shift in J — Ks color due to one magnitude of extinction at V. 
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Fig. 8. — The sdB only distribution in J— Kg (from Figure 7) fitted by the sum of two Gaussians. 
Shaded bins were excluded when calculating x 2 an d Xr- See §3.4 for further discussion, and Table 5 
for a tabulation of the fit parameters. 
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Fig. 9. — Color-magnitude diagrams of J (top) and K$ (bottom) vs. {J—K$)- The sdBs are plotted 
as circles and the sdOs as crosses ("x"). The K$ vs. (J—Ks) diagram shows the cut made to get an 
approximately volume limited sample (thick line), as well as the limits (thin lines) for companions 
with colors corresponding to the Pop I main sequence spectral types: GO, KO, K5, and MO (these 
limits are the lines from Figure 3 transformed into color-magnitude space by assuming an apparent 
magnitude for the subdwarf of Ks = 15.75). 
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Fig. 10. — Same as Figure 2, except displaying only hot subdwarfs in the "volume limited" sample; 
sdBs (176 in the left panel, 180 in the right) are denoted with circles and sdOs (62 in both panels) 
with crosses ("x"). 
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Fig. 11. — Same as Figure 7, except showing only hot subdwarfs in the "volume limited" sample. 
A tabulation of the numbers and fractions of single and composite hot subdwarfs can be found in 
Table 6. The arrow indicates the shift in J—Ks color due to one magnitude of extinction at V. 
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Fig. 12. — The "volume limited" sdB-only distribution in J — K$ (from Figure 11), fitted by the 
sum of two Gaussians. Shaded bins were excluded when calculating x 2 an d Xr- See §4.2 for further 
discussion, and Table 7 for a tabulation of the fit parameters. 
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Table 1. Updated Coordinates and Notes for All Subdwarfs with a 2MASS 2IDR Identification 

(ordered by RA). 



No. 1 


KHD ID 2 


RA (J2000) 


Dec (J2000) 


Typo 


Notes 3 




SIMBAD ID 


4 


1525 


CSHK 22957-23 


00 01 32.285 


-05 19 16.76 


sdB 






BPS CS 22957- 


23 


1526 


CSHK 22876-28 


00 01 37.625 


-35 39 53.09 


sdB 






BPS CS 22876- 


28 


1527 


PG 23595+1942 


00 02 08.520 


+ 19 59 13.04 


sd 






PG 2359+197 




1 


SB 7 


00 03 24.371 


-16 21 06.34 


sdB 


#3 (CSHK 29517- 


41) 


(BPS CS 29517- 


-41) 


6 


TON S 137 


00 04 30.992 


-24 26 21.15 


sdO 


#5 (CSHK 29503- 


-9) 


(BPS CS 29503- 


-9) 


8 


CSHK 22876-35 


00 05 56.689 


-34 53 06.49 


sdB 






BPS CS 22876- 


35 


10 


TON S 140 


00 06 46.257 


-27 20 53.60 


sdOB 


#11 (CSHK 29503 


-2) 


(BPS CS 29503- 


-2) 


12 


PHL 678 


00 07 33.785 


+ 13 35 58.70 


sdB 










19 


TON S 144 


00 10 06.925 


-26 12 56.94 


sdO 


#17 (CSHK 29503- 


-25) 


(BPS CS 29503- 


-25) 


22 


PHL 2726 


00 12 27.751 


+03 54 31.64 


sd (B) 











Note. — The complete version of this tabic is in the electronic edition of the Journal. The printed edition contains only 
a sample. 

■'■The sequence number, assigned by us, is not intended to be used as a new designator in the literature, but is useful in 
navigating the KHD catalog and cross-referencing this table with Table 2. 

2 The first identifier that appears in KHD for a given entry. 

^#NNNN (Identifier) = indicates a duplicate entry sequence number and the first identifier listed for the duplicate in 
KHD (the duplicate entry is not otherwise listed in this table); a — #87, Hebcr ct al. (2002) report that this pair is resolved 
into three components with HST: the sdO is the NW star of the "pair," while the SE "companion" star resolves into two 
components separated from the sdO by 3^37 and 4''48. b — #284, visual double, companion is M-typc (KHD); c — #293, 
Ulla ct al. (2001) report two close companions ("B" located < 2^4 E. and "C" located ~ 3^5 NW); d = #1337 and 1322 
were listed with the same name (CSHK 29495 — 21) by KHD when in fact they are two different objects: #1322 is named 
correctly, while #1337 is not — #1337 should read CSHK 29495-63 in KHD. 

4 The identifiers listed in the last column are consistent with SIMBAD nomenclature. For example, the designator 
CSHK 22876-35 used by KHD for star #8 in this table is equivalent to the SIMBAD designator BPS CS 22876-35, 
and the designator PG 23595+1942 used by KHD for star #1527 in this table is equivalent to the SIMBAD designator 
PG 2359 + 197. If no other ID is given in the last column, the identifier from KHD (given in column 2) is consistent with 
SIMBAD. In this way, one SIMBAD-compliant designator is available for cross-reference. The exceptions arc the "HS" star 
#390 HS 1000+4704, the "SC" star #1165 SC 1721-3338 from Reid ct al. (1988), and the "LBQS" stars from Berg ct al. 
(1992, #41, 565, and 1405) which arc not presently recognized by SIMBAD; although SIMBAD does use the LBQS designator 
for quasars from Hewett et al. (1995). In the case of a duplicate entry, if the duplicate's identifier listed in column 6 is not 
consistent with SIMBAD notation, then a SIMBAD-compliant name for the duplicate is given enclosed in parentheses. 



Table 2. Magnitudes and Colors for All Subdwarfs (ordered by J2000 RA). 



No. 1 


ID 2 


J 




H 


v(Hf 


K S 


<?(Ks) 3 


J~K S 3 


a (J~K s f 


V (or yf 


B- V 


Type 


Notes 5 


1525 


CSHK22957-23 


12.684 


0.023 


12.305 


0.035 


12.209 


0.030 


+0.475 


0.038 


13. 7B 




sdB 




1526 


CSHK22876-28 


15.038 


0.044 


15.063 


0.099 


15.169 


0.199 


-0.131 


0.204 


14. 3B 




sdB 




1527 


PG23595+1942 


16.238 


0.097 


16.246 


0.230 


17.103 




<-0.865 




15.44P 




sd 




1 


SB 7 


12.686 


0.032 


12.631 


0.036 


12.579 


0.034 


+0.107 


0.047 


12.67j/ 


-0.161 


sdB 


c 


6 


TONS 137 


14.534 


0.037 


14.714 


0.064 


14.765 


0.113 


-0.231 


0.119 


13.87j/ 


-0.273 


sdO 




8 


CSHK22876-35 


14.790 


0.035 


14.534 


0.060 


14.361 


0.085 


+0.429 


0.092 


15. 4B 




sdB 




10 


TONS 140 


14.543 


0.031 


14.687 


0.065 


14.863 


0.101 


-0.320 


0.106 


13.97a 


-0.253 


sdOB 




12 


PHL 678 


13.301 


0.042 


13.313 


0.040 


13.343 


0.048 


-0.042 


0.064 


13.032 


-0.076 


sdB 


s, AEA 


19 


TONS 144 


12.782 


0.030 


12.568 


0.024 


12.557 


0.030 


+0.225 


0.042 


12.8i/ 


-0.15 


sdO 




22 


PHL 2726 


13.329 


0.030 


13.397 


0.030 


13.426 


0.046 


-0.097 


0.055 


13.16a 


-0.125 


sd (B) 


s 



Note. — The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample. 

■'■The sequence number, assigned by us, is not intended to be used as a new designator in the literature, but is useful in navigating the KHD catalog and cross-referencing 
this table with Table 1. 

2 Idcntifiers listed as they appear in KHD. SIMBAD-compliant names are given in Table 1. 

3 Rcported errors arc ltr. If no error is listed, then the star was not detected in that band, and the quoted magnitude is a 95% confidence upper limit, for J — K g the 
reported value is the maximum color index when an upper limit on Kg is reported, and it is the minimum color index for an upper limit on .7. 

00 

y = Strtimgren y magnitude listed instead of Johnson V; P, J, or B — Photographic magnitude listed for V; colon (:) — uncertain magnitude. 0"5 

°s — reported single sdB or sdO star used to determine the single hot subdwarf locus (sec §3.1 and Figure 1); c — reported composite sdB or sdO star used to determine 
the composite hot subdwarf locus (see §3.1 and Figure 1); CC — 2MASS magnitudes listed are for a nearby companion only (the hot subdwarf was not detected); n — 
excluded from the sample because of a "non-subdwarf " classification; AEA — V and B — V from Allard ct al. (1994J; BBL — #1076 Sz 1095 had synthetic measurements 
of b and y only, sec §2.3; LAN = V and B — V from Landolt (1992); MEA = #320 & 1327, Stromgrcn uvby from Mooncy ct al. (2000); TYC = V and B — V for #234 are 
converted from Tycho magnitudes (sec §2.3): V T = 10.194, (B T — V T ) = +0.111 (ESA 1997); WMG = V and B-V from Williams ct al. (2001); a = #87, Hcbcr ct al. 
(2002) report this pair is resolved into three components with HST: the sdO is the NW star of the "pair," while the SE "companion" star resolves into two components 
separated from the sdO by 3' ; 37 and 4^48; b — #203, 2MASS noted a horizontal stripe due to bright star affecting the H magnitude only, however, no stripe or nearby bright 
star could be seen in the 2MASS image; d — #284, visual double, companion is M-typc (KHD) (only the companion was detected in the 2MASS 2IDR, thus it was excluded 
from the analysis); e — #293, Ulla et al. (2001) report two close companions ("B" located < 2.4 E, and "C" located ~ 3. 5 NW); / — #311, Stromgren photometry is 
from Wcscmacl ct al. (1992) but was not listed in KHD even though observations for other stars from Wcscmael ct al. (1992) were included in KHD; g — #428, IR excess 
possibly due to a wind (Thcjll ct al. 1995). 
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Table 3. Breakdown of the Hot Subdwarf sample discussed in §2. 



Sub-Sample Number of Objects 

1527 Objects in KHD 

Unrecovered 57 

Duplicate 11 

Input to 2MASS 2IDR 1459 

593 IDs from 1459 input to 2MASS 2IDR 

Subdwarfs 571 

Non-subdwarfs 22 



571 Hot Subdwarfs from 2MASS 2IDR 



With 2MASS J and K s 


385 


With 2MASS J but 95% upper confidence on K s 


180 


Detection of close companion to subdwarf only 


6 


With V (or y) 


348 


With BV (ox uvby) 


339 


385 with J and K s 


sdB/sd 


294 


sdO/sdOB 


91 


With V (or y) 


277 


With BV (or uvby) 


273 


277 with V (or y), J, and K s 


sdB/sd 


206 


sdO/sdOB 


71 


273 with BV (or uvby), J, and K s 


sdB/sd 


202 


sdO/sdOB 


71 
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Table 4. Breakdown of hot subdwarf colors based on J—Ks for the entire sample discussed in 

§3.3. 





Numbers 




Percentages 1 


Group 


Total 


Blue 2 


Red 3 


Blue 2 


Red 3 


sdB 


294 


171 


123 


58 ± 6% 


42 ± 6% 


sdO 


91 


48 


43 


53 ± 10% 


47 ± 10% 


Both 


385 


223 


162 


58 ± 5% 


42 ± 5% 



Errors reported are 95% confidence sampling errors of a bi- 
nomial distribution. 

2 "Blue" is denned as J-K s < +0.05 for the sdBs and the 
total sample, and J—Ks < —0.05 for sdOs. 

3 "Red" is denned as J - K s > +0.05 for the sdBs and the 
total sample, and J—Ks > —0.05 for sdOs. 
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Table 5. Parameters for the Gaussian fit discussed in §3.4 and shown in Figure 8. 



Parameter Blue Gaussian Red Gaussian 



Center 


-n 1 ^8 +0 - 004 

U.lOO_ 005 


, n oil +0.014 


Amplitude 


35^ 


15 ±2 


Width 


n nqi +0 - 009 

u.u»±_ 002 


n i rro+0.008 
U.lOO_ 020 



Integral of Best Fit: 278 



Errors reported are 95% marginal confidence inter- 
vals. 
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Table 6. Breakdown of hot subdwarf colors based on J—Ks for the "volume limited" sample 

discussed in 54.2. 





Numbers 


Percentages 1 


Group 


Total 


Blue 2 Red 3 


Blue 2 


Red 3 


sdB 


228 


170 58 


75 ± 6% 


25 ± 6% 


sdO 


72 


45 27 


62.5 ± 11% 


37.5 ± 11% 


Both 


300 


219 81 


73 ± 5% 


27 ± 5% 



Errors reported are 95% confidence sampling errors of a binomial 
distribution. 

o 

"Blue" is defined as J — Ks < +0.05 for the sdBs and the total 
sample, and J~Ks < —0.05 for sdOs. 

3 "Red" is defined as J-K s > +0.05 for the sdBs and the total 
sample, and J ~Ks > —0.05 for sdOs. 
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Table 7. Parameters for the Gaussian fit to the "volume limited" sample discussed in §4.2 and 

shown in Figure 12. 



Parameter 


Blue Gaussian 1 


Red Gaussian 1 


Center 

Amplitude 

Width 


-n 1 ^4+ - 013 

34 ±1 
u.uyo_ 004 


+0 27fi +0 - 002 
10 ±1 

n lf)1 +0.002 
U.1U±_q 005 


Integral of Best Fit: 


214 



Errors reported are 95% marginal confidence inter- 
vals. 



